Is dark matter a phantom ? 
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Abstract 

We show that Einstein's equations in a non-standard gauge have 
vacuum solutions with an asymptotically flat rotation curve as it is 
observed in the dark halos of galaxies. Introducing a material disk 
into this model we find a matter density in accordance with the Tully- 
Fisher relation. 
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Choosing spherical coordinates (t, r, <f>) we write the static, spherically 
symmetric line element in the form 



ds 2 = g ap dx a dx p = e a dt 2 - e b dr 2 - r 2 e c (cM 2 + sin 2 M<p 2 ). (1) 

where the metric functions a(r),b(r), c(r) depend on r only, and the speed of 
light is = 1. In (1) one usually sets c = by redefining the radial coordinate 
r. This is the standard gauge which is assumed in the proof of Birkhoff's 
theorem and leads to the Schwarzschild metric. 

Using this gauge to describe galaxies one is then forced to postulate dark 
matter in order to get the observed flat rotation curves right. However, one 
may have doubts whether this is the right idea because physically r is defined 
by the measuring process of the astronomer and with the above redefinition 
one loses contact to physics. To be on the safe side we keep c(r) ^ 0, 
but then, clearly, the field equations do not determine all three functions 
a, b, c uniquely. To choose a physical gauge fixing condition we consider the 
circular velocity V(r). From the geodesic equation on finds the following 
expression for V(r) ([2], eq.(2.14)) 

V 2 = (2) 
9tt 1 ^ 

where the T's are the Christoffel symbols for the metric (1). We now take this 
relation as our non-standard gauge fixing condition [1]. That means V(r) 
must be given, it cannot be predicted from the vacuum equations alone. 
As we shall see it is the other way around: V(r) determines the metric (1) 
uniquely. 

The non-vanishing components of the Ricci tensor for the metric (1) are 
the diagonal elements [1] 

R tt = V V + \a' 2 - \a'b' + a'c' + *a') (3) 

1 b' 4 a 12 c' 2 2 

R rr = -±(a? + 2c") + b -(a' + 2c' + ±) - ^ - *- - ±J (4) 

R m = e c - b [-l - ^c" - r(2c' + - jc'(a' -b' + 2c')} + 1 (5) 

= sin 2 t?^, (6) 
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the prime always denotes d/dr. Then the Einstein's equations without mat- 
ter can be reduced to the following three differential equations 



G t 



^a—b 



L' 2 + h'c' + -(b'-^')\+\(e' 
4 2 r J r 2 



a -c gO— b\ 



1 u , !/_/ 



„/2 



G rr = ±oV + -(a' + c') + ^- + - e b - c )= 



(7) 

(8) 



GW = ye c b 



a " +c "_I(6'-a'-2c , ) + -(a /2 -a / 6 , +aV-6V+ C /2 ) = 0. (9) 
r 2 



It is not hard to see [2] that there are only two independent field equa- 
tions. Indeed, using (8) b can be expressed by a and c. Eliminating b there 
remains one second order differential equation for a and c: 



c" = — c' + - )+-~+a'J + — + -(a' + c' . 
a' V r j r A 2 r 

Introducing the new metric function 

f(r) = c(r) + 2 log - 



(10) 



(11) 



where r c has been included for dimensional reasons, equation (10) assumes 
the simple form 



7-^ = a+ T 

This can immediately by integrated 

/' / 
log — = a + - + const. 
a' 2 



(12) 



(13) 



On the other hand the circular velocity squared (2) expressed by the metric 
functions simply becomes 



V 2 (r) = u 



r 



Using this in (13) we have 



C T 

a + - + log — = - log u. 

2 r c 

Differentiating and eliminating d by means of (11) and (14) 

2 a' 2 



f 



u 



(14) 



(15) 



(16) 
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(17) 



This gives the first diagonal element of the metric 



9u = e a = 



K a 



(18) 



1 + 2u 



where K a is a constant of integration. Then from (15) we get 



9M = 




l + 2u\2K c 
u ' r 2 ' 



(19) 



where K c is another integration constant which contains r c . Finally, exp6 
follows from (8) 



We have succeeded in expressing the metric by the circular velocity squared 
u{r). If we choose u = r s /2{r — r s ), we recover the standard gauge c = 
and the Schwarzschild metric. But now also flat rotation curves are possible 
without dark matter. As the ether is superfluous in special relativity, so 
seems to be dark matter in general relativity. From (18-20) we are able to 
predict other observable quantities which can be computed from the metric, 
for example lensing data [3]. In this way the theory can be tested. 

Another test is to introduce normal matter located in the equatorial 
plane which simulates a spiral galaxy. There exists a simple method to 
construct solutions of Einstein's equations with such a disk. This is the so- 
called displace, cut, and reflect method which goes back to Kuzmin [4] and 
since then was used and modified by many authors (see [5] and references 
given there). If the equatorial plane is z = in cylindrical coordinates, the 
matter density is singular ~ S(z) and the metric g^ u has jumps in its first 
derivatives. The mathematical basis in this situation is Taub's theory of 
distribution valued curvature tensors [6]. The Einstein tensor then has a 
singular part which is related to the distribution valued energy-momentum 
tensor by 




(20) 




= Kt^S(z), 



(21) 



where 



R = g' 




8ttG 



4 



If the jumps h^u of the normal derivatives of g^ u are defined by 

[W] - |^ , -fe = n^, (22) 



<9x° 



where + and — mean the limiting values from both sides of the singular 
plane with normal vector n a , then t^ u can be calculated from (6], eq.(6-2))): 



-H, = « 2 ((^^)(i-^)- 

-(9r-^)(r-™))bor, (23) 

where n 2 = n a n a . 

We go over to cylindrical coordinates (t, R, z, 4>) 

r 2 = R 2 + z 2 , z = rcosi3, sin$ = — . 

r 

Then the metric (1) assumes the following non-diagonal form 

ds 2 = g iiv dx p 'dx v 

with 

p2 r 2 p2 _2 

_ n i z n z i 
9oo = e , 5n = o~e n e ' 522 = n~e ^e" 

512=521 = -2^ (e fe -e c ), g 33 = -R 2 e c . (24) 

For simplicity we still use r, but our admissible coordinates are x 1 = R,x 2 = 
z. 

We now apply the displace, cut, and reflect method following the pro- 
cedure of Voigt and Letelier [5]. We take the metric (24) with a,b,c given 
by (18-20) in the half space z > d > 0, displace it to z = and reflect it 
for z < 0. This produces the finite jumps in the z-derivatives of g^ v . The 
whole procedure is equivalent to the transformation z — > The normal 

vector is n M = (0, 0, 1, 0) = 5^ and 

n v = <f»n ll =!? 2 , n"n u = g 22 . 

The jumps (22) in the normal derivatives on z = which we need are equal 
to 

foil = [511,2] =0n—- -J e (25) 
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t r , , 2d 

033 — [533,2] — #33 ~> 

where the prime always means d/dr keeping z and R constant. Now from 
(23) we find the energy density 

(26) 

with D = ^11522 - 9i2- Using 

t 2d/ , 2 c x 2d , 

hi = — [9ii e , 633 = — 533c, 

we finally obtain 

Here we have to put z = everywhere which gives r 2 = i? 2 + d 2 . 

Now we must specify the circular velocity squared u(r) in order to fix 
the metric. We are particularly interested in the case of an asymptotically 
flat circular velocity which in the usual terminology corresponds to a dark 
halo. Therefore we assume u{r) of the form 

u(r) = u^ t + ^ + 0(r~ 2 ) (28) 

for large r. Then it follows from (18-20) 

e? = K a + 0{r- 1 ), e b = ^J+0(r- 5 ) (29) 



L. 

r 

where by (19) 

" x flat = ^flat 

Using this in (27) the leading order comes from the last term 



e C = -3 + °(^ 3 ) ( 3 °) 



Lc ~ < = V a ~t- (31) 



2r7 r 2 

*8 = ^^(l + 0(/2- 1 )). (32) 

This is proportional to the density of normal matter because we consider a 
static energy-momentum tensor. Taking (31) into account we find that 

4-<t-V R \ t (R) (33) 
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for large R. This is in accordance with the baryonic Tully-Fisher relation for 
galaxies [7] [8], which states that the total baryonic mass M is proportional 
to Vg 4 at . In fact, one can show [9] that the contribution of the inner part 
R < Ri of the disk can be made arbitrarily small compared to the outer 
part between R\ < R < i? 2j say. We emphasize that M is obtained from 
i[j by integrating with the Euclidean surface measure R dR d(j), because this 
is what astronomers are doing when they determine M from luminosity 
measurements. 

The radial pressure t T r vanishes because G r r (8) does not contain a second 
derivative. Therefore our model must be interpreted as a dust disk with 
purely azimuthal stresses. This is not very realistic and it remains to be 
investigated whether the Tully-Fisher relation is a generic property for more 
physical galaxy models. But at present the answer to the question in the 
title is a careful "may-be yes" . 
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